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Leptogenesis constitues a very simple scenario to achieve the baryon asymmetry that we observe today. It 
requires only the presence of right handed neutrinos (which arise very naturally in many extensions of the Standard 
Model) and depends crucially on the mechanism responsible for their production. In particular, when their mass 
exceeds the inflaton mass (10^'^ GeV in chaotic inflation) only non perturbative production can occur. It is shown 
here that non perturbative production of fermions in an expanding Universe is a very efficient mechanism up to 
masses of order of flO^^ — 10^**) GeV, and that it can be easily applied to solve the baryon asymmetry problem. 



1. Introduction 

Inflation is at present the most natural and 
accepted solution to many cosmological puz- 
zles, such as the flatness, the horizon, and the 
monopole problems. However, the Universe af- 
ter inflation looks very different from the one we 
observe today, since the huge expansion leaves it 
empty of radiation and matter. 

The only form of energy which is not diluted 
away after this phase is the vacuum energy of the 
field (or fields) which leads inflation itself, namely 
the inflaton field. The process of converting this 
energy into the radiation and matter which we see 
around has been named reheating and it has been 
the object of intense studies in the past decades. 

Apart from this energy conversion, the pro- 
cesses which govern the first instant of the Uni- 
verse must also favor the production of matter 
over antimatter. Indeed, on one hand direct ob- 
servations put strong limits to the presence of an- 
timatter in our cluster, while on the other consid- 
erations about how the light elements abundances 
were formed when the Universe was about 1 MeV 
hot lead us to conclude that the difference be- 
tween the number density of baryons and that of 
antibaryons is about 10^ ^^ if normalized to the 



entropy density of the Universe. 

Until now, many mechanisms for the genera- 
tion of the baryon (B) asymmetry have been pro- 
posed (for a recent review, see [||). 

Grand Unified Theories (GUTs) unify the 
strong and the electroweak interactions and pre- 
dict baryon number violation at the tree level. 
They are — therefore — perfect candidates for 
a theory of baryogenesis. There, the out-of- 
equilibrium decay of superheavy particles can ex- 
plain the observed baryon asymmetry. 

Another plausible scenario is the theory of elec- 
troweak baryogenesis (in the context of the Min- 
imal Supersymmetric Standard Model), where 
baryon number violations take place at the quan- 
tum level due to the usoppressed and baryon 
number violating sphalerons in the hot plasma. 

Sphaleron transitions leave unchanged the 
combination B ~ L (where L is the lepton num- 
ber) and play the role of converting an initial lep- 
ton asymmetry to a final mixture of nonvanish- 
ing B and L . This is the key idea of the theories 
of baryogenesis via leptogenesis. Once a lepton 
asymmetry is produced, sphaleronic transitions 
convert a fraction of it into baryon number. In 
the Standard Model the following relation for the 
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final amount of the baryon asymmetry liolds: 
8nf + Ann 



B = 



(B-L) 



(1) 



22n/ + rSuH , 

wtiere nn is tlie number of Higgs doublets and Uf 
is the number of fermion generations. 

Adding right-handed Majorana neutrinos to 
the SM breaks B — L , and the primordial lep- 
ton asymmetry may be generated by their out- 
of-equilibrium decay. This simple extension of 
the SM can be embedded into GUTs with gauge 
groups containing SO{10). Heavy right-handed 
Majorana neutrinos are particularly welcome to 
explain the smallness of the light neutrino masses 
via the see-saw mechanism |3|. In most of 
the models proposed so far to generate neutrino 
masses, the latter are of the Majorana type, which 
implies the existence of interactions which violate 
the lepton number and thus render leptogenesis 
particularly attractive. 

The leptogenesis scenario depends crucially on 
the details of the reheating process after inflation 
and on the production mechanism that was re- 
sponsible for populating the Universe with right- 
handed neutrinos. 

The simplest way to envision this reheating 
process is if the comoving energy density in the 
zero mode of the inflaton decays perturbatively 
into ordinary particles, which then scatter and 
thermalize to form a thermal background Ppl. 
Of particular interest is a quantity known as the 
reheat temperature, which is given by 

Trh = v/r^Mp , (2) 

where F^ is the inflaton decay rate. 

In order to produce a sufficient amount of 
heavy right-handed neutrinos with mass Mj^ , it 
is commonly assumed that Tr// ^ M^ . 

However, T^h is constrained by the so called 
gravitino problem m to be less than about 
IQi" GeV. 

This bound increases if the heavy neutrinos are 
produced directly through the inflaton decay pro- 
cess. This would be kinematically accessible only 
if0 

Mn < M^ - 10" GeV . (3) 

^The last equality is required by the observed density and 
temperature fluctuations. 



The outlook for leptogenesis might be bright- 
ened even further with the realization that re- 
heating may differ significantly from the simple 
picture described above |0,|[ . In the first stage of 
reheating (namely in the first dozen or so oscil- 
lations), called preheating [|7|], nonlinear quantum 
effects may lead to extremely effective dissipative 
dynamics and explosive particle production, even 
when single particle decay is kinematically forbid- 
den. 

It has been already shown that preheating 
plays an extremely important role in baryogene- 
sis, in particular for the production of superheavy 
GUT bosons PJlC|]. The application of preheat- 
ing that will be instead presented here is a very 
new one M and it consists in enriching the lepto- 
genesis scenario with the possibility of producing 
right-handed neutrinos which are too heavy to be 
generated in the usual perturbative ways. 

In what follows it will be shown (both with 
numerical calculations and with some analytical 
estimates |^ that the production of a generic su- 
perheavy fermionic field during preheating is a 
very efficient mechanism up to masses of order of 
10^* GeV. As a consequence, it will be pointed out 
that a sufficient L asymmetry can be generated if 
the fermionic field is a right-handed neutrino. 

2. Production of fermions during preheat- 
ing 

The coupling inflaton — fermions gives to the 
latter the effective time dependent mass 



TO(i) = mx + g(f>it) 



(4) 



where mx is the bare fermionic mass. 

This changing mass induces the creation of X 
particles. To see this, one has to diagonalize the 
hamiltonian associated to the Dirac equation 



a 



.Ih-,' 



m] X ^0 



(5) 



that is written here for an expanding flat Universe 
{H is the Hubble constant) and in conformal time 

V ■ 



^New analytical expressions with respect to the paper llj 
— where the reader is addressed for a more general treat- 
ment of leptogenesis — will be presented here. 



Once the diagonalization is performed (see [y 
for details) the number of produced fermions is 
given by 



n {rj) : 



T^^a^iv) Jo 



dkk Uf. (77) 



"fc iv) - 



2lu 



-(6) 



2 CO — 2klle (u*^u^) ~ am (1 ~ 2u*^u+) 

where uj'^ = k'^ + m (t) a^ , fc is the momentum 
of the produced particle, and a is the scale fac- 
tor of the Universe. The functions u± are, re- 
spectively, the positive and negative eigenfunc- 
tions of 7o which enter in the decomposition of 
the fermionic field. They satisfy the uncoupled 
equations 



— :t + uj zLi(a'm + am') 
drj'^ 



u±{k) = 0, 



(7) 



with initial conditions at time rj — given by n 



..(0, - ^1(1.!^) 

m'-i-(O) = iku^(0) ^ iamu±(0) 



(8) 



A measure of the strength of the coupling in- 
flaton — fermions is given by the parameter 



q = g^<j>^iO)/Aml 



(9) 



where </> (0) ~ 0.28Mp is the value of the inflaton 
at the end of the inflation and m^ ~ lO^'^ GeV is 
the inflaton mass. 

In Fig.Q^ the results of a numerical integration 
of the equations given above are reported. For 
a wide range of mx and q the value of the ratio 
px/p is shown, where px = rnx-nx is the energy 
density transferred to the fermions and p is the 
initial energy density stored in the inflaton field. 

One can notice the presence of the cut-off 
"Jmax ~ <7^'^/2 in the fermionic mass. For mx 
higher than this value the production is substan- 
tially zero, while for lower values it increases lin- 
early with q and is independent from mx ■ This 
last statement seems not to be completely con- 
firmed by FigJ^, where the value px increases 

*Note that with these values n (0) = . 
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Figure 1. Ratio between the energy density of 
the produced fermions and of the inflaton field, 
for many values of the fermion mass and of q . 



slowly with mx ■ However, as the analytical esti- 
mate below will confirm, the scaling stated here 
is exact, while the slope that one sees in FigJ^ for 
high q and small mx is a numerical artifact. The 
problem is that — due to limited computational 
resources — px reported in Fig.|l| is evaluated af- 
ter 20 oscillations of the inflaton field. However, 
the number of oscillations of <j) that one needs for 
px to saturate increases with q , and a slope in px 
appears for those values for which the saturation 
has not yet taken place. 

To go on, it is important to understand the 
behavior of Uk , that is the number density of 
fermions produced with a given momentum k . 
All these spectra exhibit the qualitative behavior 
that is reported in Fig.|| for m,x = 100 m,^ and 
for two values of q . 

Also here one can find a threshold fcniax in 
the production. For k < fcmax , nk approxi- 
mately saturates to 1 (filling the whole allowed 
Fermi sphere), while for momenta bigger than 
fcmax the production does not occur. Thus one 
can estimate [see Eq. (ph] 



Px = mx ■ nx oc mx 



(10) 



As it is shown in Fig.0, fcmax increases with 
q . More precisely, the analytical estimate that 
will be given below indicates fcmax cc (q/mx) 



rewrites 




Figure 2. Spectrum in momenta of the produced 
particles for a given fermion mass and two differ- 
ent values of q . 



From this value and from Eq. (nm one immedi- 
ately gets the scaling px oc q stated above. 

To understand the scaling of fcmax, one can 
differentiate Eq. (ph: 



drj 



duj (P 
drj 



Vdr/2 



(771a) 



(11) 



This expression suggests that n^ changes signif- 
icantly when the total mass m approaches zero 
and when the variation of the frequency uo is non- 
negligible. Indeed, our numerical results indicate 
that the production of fermions occurs for very 
short intervals about the zeros of m . The esti- 
mate that is here presented consists in neglecting 
the last factor in Eq. (^l]) and, as it is custom- 
ary done for preheating of bosons, in assuming 
that the production occurs only for Iw'/w^l > 1 , 
that is in a regime of non adiabatic change of the 
frequency R. 

In terms of the number N of oscillations of 
the inflaton field, the non adiabaticity condition 



^The results that will be presented here can be obtained 
(in a more lengthy way) also without this assumptions and 
hence do not depend on the approximations made here. 



J \ 2/3 
m dm \ 



„ < I — m" 

-'2 I 271 dN ) 



(12) 



where masses and momenta are given in units of 
the inflaton mass. 

The cut-off in the production for mx > Wmax 
can be now understood by rewriting Eq. (0) in 
terms of N : 



m = mx + '2q 



1/2. 



0(0) 



mx 



' ttN 



cos2 7riV.(13) 



At large mx the first zero of m can be potentially 
reached at A^ = 1/2. However, for m.x greater 
than 2g-'-/^/7r, the total mass m, never vanishes 
and Eq. ( [l2[) is never satisfied Q. 

For mx < ?Timax , the total mass vanishes at 
the points Nq given by (iVmax is the greatest Nq ) 

„l/2 „l/2 

No = -- cos27rAro , N^ax ^ . (14) 



TT mx 



IT mx 



To find fcniax , one has to maximize r.h.s. of 
Eq. (p^. The maximum occurs when the two 
terms in r.h.s. are of the same order. Thus, a 
good estimate for fcniax can be achieved by solv- 
ing 



m,m 

2-n 



2/3 



(15) 



expanding m and m! for N close to Aq . Doing 
so, one gets 



k-^am^ |cos27rAro|^/^9^/^A^o^^. 



(16) 



Since k grows with A'o , k^ax has to be eval- 
uated in Nmax ■ Thus the particles with the 
highest momentum are generated in the last zero 
Nmax of TO . This gives the stated result 






1/3 



(17) 



®While numerical coefficient based on the second ap- 
proximate equality in the r.h.s of Eq. (hsl) should not 
being taken seriously, numerical integration shows that 
2(/)/0(O) ~ —0.5 in the point of first minimum of </>, which 



gives the cut-off in the proper place, mx 



u/2 



/2. 



3. Baryogenesis 

The previous section shows that preheating can 
produce fermions with the energy density 

-.1/2 



PX 
P<t> 



IQ- 



q up to AIx 



q 



(18) 



This resuh can be useful for leptogenesis if one 
identifies the X particles with right-handed neu- 
trinos N and if these heavy neutrinos decay in a 
CP violating manner into left handed neutrinos 
plus higgs particles. 

Eq. (18) is obtained with the tacitly assump- 
tion that the produced fermions are stable. How- 
ever this estimate does not apply if the right 
handed neutrinos have a decay lifetime larger 
than the typical time-scale of the inflaton oscilla- 
tion TTlT . 

Another requirement that must be put is that 
the right handed neutrinos do not annihilate into 
inflaton quanta more rapidly than they decay in 
left handed neutrinos, otherwise the produced 
lepton asymmetry would be negligible. 

Both the above requirements can be typically 
satisfied V\ for masses M^ < 10^^ GeV and for 
couphngs q < 10^° . 

Once the right handed neutrinos decay and 
once the lepton asymmetry produced by the 
decay is converted in baryon asymmetry by 
sphaleron transitions, Eq. ( [l8| ) gives ||l| 



3^"^ = 



3-10' 



Tr. 



10i"GeV; V M, 



/IQi^GeVN / q 



iN 



) VioioJ 



,(19) 



where T^h is the reheating temperature and e is 
the amount of CP violation in the decay of the 
heavy neutrinos n. 

Eq (|l9|) reproduces for very natural values of 
the parameters the observed baryon asymmetry 

'^One must also ask that lepton number violating scatter- 
ings do not deplete the lepton asymmetry generated bye 
the d eca y of the right handed neutrinos. See ||l| for details. 
*Eq (hfl) holds if the energy density of the right handed 
neutrinos never dominates over the thermal one. However, 
if the neutrinos have a very late decay they can eventu- 
ally dominate over the thermal bath, since their energy 
does not redshift with the expansion of the Universe. If 
it is the case, the decay of the neutrinos gives both the 
baryon number and the entropy of the thermal bath, and 



of the Universe. This result leads thus to conclude 
that baryogenesis via leptogenesis is a viable op- 
tion also if one considers right-handed neutrinos 
with mass higher than the inflaton one. 
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the resulting baryon asymmetry is given by 

B = 8-10'^ { ^ e , (20) 

where m is the mass of the light neutrinos. 



